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1.  INTRODUCTION 


Lightning  is  of  interest  to  the  Air  Force  for  a  number  of  reasons.  It 
is  a  threat  to  aircraft,  both  in  flight  and  on  the  ground.  It  is  a  threat  to 
ground  installations,  particularly  antennae  and  associated  electronic  equip¬ 
ment  and  power  transmission  lines.  It  is  a  potential  spurious  trigger  of  the 
Nuclear  Detonation  Detection  System. 

In  addition  to  the  reasons  given  above  for  the  importance  of  natural 
lightning  in  its  own  right,  an  understanding  of  natural  lightning  is  necessary 
as  a  basis  for  understanding  nuclear  lightning. 

In  order  to  characterize  the  threat  posed  by  lightning,  the  important 
parameters  are  the  (peak)  current  and  the  rate  of  current  rise.  Unfortu¬ 
nately,  the  observations  of  lightning  to  date  do  not  provide  a  basis  for  con¬ 
fidently  predicting  or  bounding  these  parameters.  To  rectify  this  situation, 
in  an  ongoing  program  involving  the  Air  Force  Weapons  Laboratory  (AFWL)  and 
New  Mexico  Institute  of  Mining  and  Technology  (New  Mexico  Tech)  Langmuir 
Laboratory  facility  on  South  Baldy  Peak,  New  Mexico,  Mission  Research  Corpora¬ 
tion  (MRC)  has  begun  a  program  to  provide  simultaneous  observations  of  cur¬ 
rent,  brightness,  and  radiated  fields  of  a  lightning  stroke.  The  results  of 
these  initial  efforts  are  described  in  this  report.  These  results  are  unique, 
and  provide  valuable  new  insights  into  natural  lightning.  For  the  first  time, 
there  is  evidence  that  the  return  stroke  velocity  is  not  a  constant  (typically 
reported  in  the  literature  as  one-third  that  of  light  with  a  maximum  observed 
value  of  about  half  the  velocity  of  light),  but  is  initially  almost  equal  to 
the  velocity  of  light  and  decreases  with  propagation  up  the  channel.  In  addi¬ 
tion,  for  the  first  time,  there  are  correlated  brightness  and  current  observa¬ 
tions.  These  are  in  good  (but  not  perfect)  agreement  with  computer  calcula¬ 
tions  of  the  optical  emissions  based  upon  the  observed  current. 

Sections  2  and  3  of  this  report  describe  the  observation  of  visible 
lignt  emission.  Computational  modeling  is  compared  to  observations  in  Section 
3.  Section  4  derives  the  relationship  between  observed  intensity  and 
intrinsic  channel  brightness  used  to  analyze  these  observations. 


Section  5  describes  the  analysis  of  the  electromagnetic  observations  obtained 
by  B-dot  and  D-dot  sensors  located  near  the  optical  sensors.  It  was  these 
observations  that  revealed  the  high  initial  return  stroke  velocities. 

Finally,  Section  6  discusses  theoretical  analysis  of  experiments  done 
at  the  Langmuir  Laboratory  on  the  use  of  Primacord  in  the  triggering  of  light¬ 
ning.  The  aim  of  the  experiment  was  to  form  a  conducting  diwswiel  rapidly  by 
detonating  a  length  of  Primacord  suspended  from  a  balloon.  Et  was  hoped  that 
this  would  provide  a  more  realistic  simulation  of  a  natural  lightning  strike 
to  a  tower  than  towing  a  conducting  wire  up  with  a  rocket.  This  method  did 
not  appear  to  provide  a  sufficiently  conducting  channel  for  an  adequate  time 
to  initiate  a  strike. 

This  work  should  serve  as  the  foundation  for  a  quantum  leap  in  the 
understanding  of  natural  lightning.  It  should  be  possiDle,  lay  applying  and 
extending  the  results  of  this  study,  to  infer  the  current  in  a  lightning 
stroke  by  photographic  observations;  this  should  be  of  value  in  understanding 
nuclear  lightning.  This  work  should  lead  to  a  better  understanding  of  return 
stroke  initiation  conditions  and  the  development  of  the  return  stroke,  leading 
to  better  prediction  of  currents  and  rates  of  current  rise.  The  threats  posed 
by  lightning  will  then  be  better  characteri zed . 


2.  DESCRIPTION  OF  LIGHTNING  BRIGHTNESS  EXPERIMENT 

2.1  INTRODUCTION 

The  aim  of  the  experiment  was  to  record  the  light  output  of  a  small 
length  of  a  lightning  channel  at  the  same  time  the  current,  electric  field, 
and  magnetic  field  were  measured.  This  would  be  the  first  time  such  simul¬ 
taneous  observations  were  made.  The  time  dependence  of  the  light  output 
was  determined  with  a  10  ns  time  resolution,  which  was  the  resolution  of 
the  Biomation  and  Transiac  digitizers  that  were  available  in  Kiva  1.  The 
photodiodes  had  a  5  ns  resolution. 

2.2  DESCRIPTION  OF  OBSERVATION  SITE 

The  observations  were  performed  at  a  facility  at  South  Baldy  Peak 
near  the  Langmuir  Laboratory  of  New  Mexico  Tech.  The  facility  includes  a 
grounded  aluminum  mast  penetrating  the  roof  of  an  underground  metal 
building,  Kiva  2.  The  mast,  which  is  20  m  long,  is  attached  at  its  base  to 
a  lightning  current  measuring  sensor  array  designed  by  EG&G.  The  array 
measures  the  electrical  current,  and  its  derivative  dl/dt  (Ref.  1).  Two 
silicon  photodiodes  used  for  measuring  optical  radiation  from  lightning 
channels  were  located  at  a  second  underground  ouilding,  Kiva  1. 
Electromagnetic  sensors  used  to  measure  two  horizontal  components  of  aB/dt 
(3  is  the  magnetic  induction)  and  the  vertical  component  of  dD/'dt  (D  is  t.ne 
electric  displacement)  are  mounted  to  the  roof  of  Kiva  1  (Ref.  2).  The 
"ground  plane"  in  the  vicinity  of  Kiva  1  and  Kiva  2  and  between  them  is 
covered  witn  a  grounded  wire  mesh.  A  map  of  the  facility  is  shown  in 
Figure  1  taken  from  Ref.  3.  The  brightness  experiment  was  designed  by 
William  Rison  of  New  Mexico  Tech  and  Carl  E.  Baum  of  the  Air  Force  Weapons 
Laboratory  in  collaboration  with  the  authors  of  this  report. 

The  geometry  for  the  brightness  measurement  is  shown  in  Figure  2.  ~'ne 
detector  views  an  area  centered  13  m  above  the  mast  of  the  current  sensor 
array  in  order  to  decrease  the  likelihood  that  metal  vapor  would  affect  the 
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light  output  from  the  channel.  Assuming  a  vertical  channel,  the  detectors 
measured  the  brightness  of  a  section  of  the  channel  that  was  12.0  m  long 
between  the  half  sensitivity  points,  9.0  m  long  between  full  sensi  tivity 
points,  and  15.1  m  long  between  zero  sensitivity  points.  The  section  of 
channel  was  centered  at  18  m  above  the  top  of  the  mast  of  the  current  sen¬ 
sor  array  mounted  on  Kiva  2.  Collimators  in  front  of  the  Lite  Mike  detec¬ 
tors  were  designed  so  that  the  horizontal  field  of  view  was  17.4  m  at  a 
distance  of  60.8  m,  twice  as  large  as  the  vertical  field  of  view;  see  Figure  3. 
The  Lite  Mike  detectors  were  configured  in  such  a  way  that  the  center  of 
the  field  of  view  of  each  of  them  coincided  with  the  center  of  the  field 
of  view  of  a  video  camera  that  was  mounted  with  them  as  a  unit.  The  video 
camera  was  used  to  aim  the  Lite  Mike  detectors.  The  metal  box  containing 
the  Lite  Mike  heads  and  the  video  camera  was  mounted  in  the  crawl  space 
adjacent  to  Kiva  1,  as  is  shown  in  Figure  4. 

2.3  RELATION  BETWEEN  LITE  MIKE  VOLTAGE  AND  OPTICAL  POWER 
PER  UNIT  LENGTH  RADIATED  8Y  THE  LIGHTNING  CHANNEL 

The  switch  on  each  Lite  Mike  detector  head  was  set  to  the  white  XI 
position,  in  order  to  match  the  impedance  of  the  coaxial  cables  and  the 
digitizers.  In  this  position,  with  the  50-fl  cables,  and  a  50-i2  digitizer, 
the  calioration  constant  for  either  detector  was  0.5  K^,  wnere  K;  is  the 
quoted  calibration  constant  in  volts  per  watt.  The  optical  power  values 
were  obtained  as  follows.  The  length  of  channel  observed,  assuming  that  it 
was  vertical,  was  i  -  12.0  m.  The  equation  for  the  flux  P  in  watts  through 
the  surface  of  the  photodiode  at  a  given  wavelength  is 

p  (X)  =  MilA  (1) 

4wR2 

and  the  voltage  is 

<  i 

V  *  —  P  (2) 

where  A  is  the  area  of  the  photodiode,  W(x)  is  the  power  radiated  per  unit 
wavelength  interval  by  a  1  m  length  of  the  channel,  P  is  the  convolution  of 
P(X)  with  the  frequency  response  of  the  photodiode,  and  the  factor  of  1/2 
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Figure  3.  C.ollimator  geometry  for  Lite  Mike  detector. 
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Figure  4.  Lite  Mike  heads  and  video  camera  were  located  in  the  crawl 
space  adjacent  to  Kiva  1  to  protect  them  from  the  weather. 
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arises  because  the  lite  mike  output  is  impedance  matched  rather  than  con¬ 
nected  to  an  infinite  impedance.  The  value  of  K \  corresponding  to  the  Lite 
Mike  with  product  serial  no.  0272  is  30.20  V/W,  and  the  photodiode  has 
area,  A  =  0.0507  cm2.  The  product  KjA  is  1.24  V*cm2/W  for  the  other  Lite 
Mike,  serial  no.  0420,  KiA  *  1.53  V«cm2/W. 


3.  CHANNEL  BRIGHTNESS  CALCULATIONS 


3.1  INTRODUCTION 

The  CHARTS  one-dimensional  magnetohydrodynamic  (MHD)  code  was  used  to 
simulate  the  optical  radiation  emitted  by  the  channel.  CHARTB  was  developed 
at  Sandia  National  Laboratories.  It  is  a  state-of-the-art  program,  con¬ 
taining  approximately  12000  lines  of  FORTRAN  code.  It  includes  a  tabular 
equation-of-state  for  air.  A  fully  consistent  calculation  of  the  current 
profile  in  the  channel  (incorporating  magnetic  diffusion  and  advection)  is 
performed  and  the  Lorentz  force  is  included.  In  air  and  other  gases,  CHARTB 
has  been  used  to  model  the^action  an  explosive  opening  switch,  to  simulate 
exploding  wires  and  to  solve  various  other  problems  with  good  results. 

CHARTB  has  been  used  to  simulate  the  evolution  of  plasma  channels  for  parti¬ 
cle  beam  transport. 

In  each  case,  the  lightning  calculations  were  begun  with  an  initial 
channel  temperature  of  1  eV  and  a  channel  with  a  1  mm  radius.  The  current 
was  specified  by  the  equation 

I  =  IQ[exp  (-t/tf)  -  exp(-t/tr)]  (3) 

where  t  is  the  time  and  the  constants  tr,  tf,  and  I0  were  selected  to 
match  Equation  3  as  closely  as  possible  to  the  measured  current. 

3.2  RADIATIVE  TRANSPORT 

CHARTB  contains  a  multigroup  radiative  transport  algorithm  using  the 
diffusion  approximation,  which  was  used  to  calculate  the  power  of  the  opti¬ 
cal  radiation  in  the  band  of  sensitivity  of  the  EG&G  Lite  Mike  detector  that 
was  used  to  measure  the  optical  power  radiated  per  unit  length  by  the  light¬ 
ning  channels.  The  radiative  transport  algorithm  used  is  described  in 
detail  in  Ref.  4.  The  spectrum  of  optical  radiation  was  divided  into  fre¬ 
quency  intervals,  or  groups.  The  group  boundaries  were  selected  to 
minimize  the  opacity  variation  within  each  group.  A  single  value  of  the 


opacity  was  assigned  to  each  group  using  the  Rosseland  mean.  This  type  of 
average  is  more  accurate  for  the  optically  thick  case,  in  which  the 
radiation  is  absorbed  and  reradiated  a  number  of  times  in  traversing  the 
calculational  grid,  than  for  the  optically  thin  case.  In  the  limit  in  which 
groups  are  narrow  enough  to  resolve  the  opacity  versus  frequency  function, 
the  algorithm  used  is  valid  in  both  the  optically  thick  and  optically  thin 
limits.  Optical  radiation  is  the  dominant  mechanism  of  energy  transfer  from 
the  hot  conducting  gas  to  the  cooler  gas  in  the  outer  part  of  the  channel. 

The  diffusion  approximation  is  also  known  as  the  method,  which 
refers  to  the  first-order  expansion  of  the  intensity  in  spherical  harmon¬ 
ics.  With  the  Marshak  boundary  condition,  which  we  used-,  this  approximation 
gives  radiant  energy  densities  that  are  accurate  in  the  optically  thick  lim¬ 
it,  but  the  energy  densities  are  overestimated  in  the  optically  thin  limit. 
However,  the  power  radiated  by  the  hot  part  of  the  channel  is  almost  inde¬ 
pendent  of  the  radiant  energy  density  for  optically  thin  bands,  and  the 
energy  in  thin  bands  that  is  reabsorbed  in  the  cooler  gas  is  a  negligible 
fraction  of  the  energy  radiated  in  the  hot  region.  Therefore,  large  errors 
would  result  if  the  calculated  output  flux  at  the  edge  of  the  grid  had  been 
used  to  obtain  the  radiated  power  in  the  band  of  interest,  but  by  summing 
the  net  radiated  power  over  all  grid  zones,  accurate  results  were  obtained. 

The  radiant  field  is  complete’y  specified  by  the  intensity  function, 

I(t,r,a,v),  where  t  is  the  time,  r  is  the  location,  a  is  a  unit  vector  in 

the  direction  of  travel,  and  v  is  the  frequency.  At  a  point  interior  to  the 

lightning  cnannel  ,  the  light  that  is  radiated  by  the  gas  in  a  differential 

volume  will  travel  away  in  all  directions.  Light  also  arrives  at  the  point 

from  all  directions,  and  any  ray  that  passes  through  the  differential  volume 

is  partially  absorbed.  The  intensity  I(t ,r , a, v)dfldv  (W/rtr$rH  )  is  the  flux  with 

frequency  between  v  and  v  +  dv  within  solid  angle  da  of  a  light  ray  with 
♦ 

direction  a  (see  Figure  z) .  The  light  ray  oasses  through  point  r. 


Figure  5.  Illustration  of  the  vectors  and  differential  elements  that 
specify  the  intensity  function. 


The  equation  of  radiative  transfer  is 


1 


+  Cj  •  7)  I  »  Sources  -  Sinks. 


The  convective  derivative  of  the  intensity  'unction  is  the  intensity  created 
minus  the  intensity  adsorbed.  The  sinks  are 

Sinks  *  c(kab$  +  kscat)  I  .  i*} 

The  intensity  loss  due  to  scattering  from  one  frequency  to  another  is 
included  in  kscat.  The  frequency-  and  temperature-dependent  absorption 
coefficient  per  unit  distance  is  kaQS.  The  source  terms  are 


SSCat  *s  the  ^tensity  gain  due  to  scattering  from  one  frequency  to 
another  and  ckatJS  B(T,v)  is  the  rate  of  increase  of  the  intensity  due  to 
the  emitted  radiation.  B(T,v)  is  the  Planck  function  for  the  blackbody 
intensity  spectrum.  The  scattering  terms  are  small  compared  to  the  absorp- 
tivities  and  may  therefore  be  neglected. 

Ignoring  the  scattering  terms.  Equation  4  becomes 

(1^.  •  7)  I  .  ckabs  (B<T,v)  -  I)  .  (7) 

In  the  optically  thick  steady-state  case  in  which  absorption  exactly  cancels 
emission,  I  *  B(T,v).  In  the  optically  thin  limit,  I  «  B(T,v)  and  the  emis¬ 
sion  rate  is  simply 

I  *  ckabs  B(T’v)  •  (8) 

This  equation  justifies  our  lack  of  concern  for  the  accuracy  of  I  for  opti¬ 
cally  thin  bands,  for  which  I  «  B(T,v).  Equation  7  was  used  to  obtain  the 
radiated  power  in  each  volume  element  rather  than  obtaining  the  radiated 
output  from  the  outward  component  of  I  at  the  edge  of  the  grid. 


The  diffusion  approximation  of  Equation  7  is  obtained  in  terms  of  the 
energy  density  of  the  radiation  field,  U  =  1/c  j  Id  3.  By  assuming  that  I  is 
nearly  isotropic,  the  diffusion  approximation  is  obtained  (Ref.  4). 

(r  It  • 7 '  Jirat)s  7)  u  •  lcabs(r“  - u).  <*) 

At  the  center  of  the  coordinate  system,  the  boundary  condition 
3U/3r  =  0  was  applied.  The  boundary  condition  applied  at  the  edge  of  the 
grid  is  a  Marshack  condition  of  the  form 


4,s^*U  =  ° 


(10) 


A  standard  tridiagonal  finite-difference  algorithm  was  used  to  solve 


Spectral  lines  were  not  treated  individually,  but  were  included  in 
the  continuum  absorpti vities.  Their  contributions  to  the  opacities  were 
therefore  somewhat  overemphasized  (Refs.  5,  6).  The  opacity  tables  of 
Ref.  6  were  used  in  the  calculations. 

3.3  CALCULATED  PROPERTIES 

The  first  lightning  strike  that  was  triggered  at  the  South  Baldy 
facility,  designated  #85220,  had  a  peak  current  of  about  17  kA  with  a  rise 
time  tr  =  2.4  x  10*7  s.  I0  and  tr  were  given  these  values  for  the 
numerical  simulation,  and  the  decay  time  was  given  a  value 
tf  =  2.5  x  10~5  s.  The  time  duration  used  for  the  simulation  was 
5  us.  The  measured  current  is  shown  in  Fig.  6  and  the  current  pulse  shape 
assumed  for  the  calculation  is  shown  in  Fig.  7. 

Radial  current  profiles,  temperature  profiles,  density  profiles,  and 
pressure  profiles  were  generated.  The  temperature  profiles  are  shown  at 
times  of  1.0  ys  ,  2.0  us  and  5.0  us  in  Fig.  8.  It  is  seen  that  the  tempera¬ 
ture  profile  has  a  dip  in  the  center  at  very  short  times.  This  shape  is 
probably  caused  by  the  assumed  initial  channel  temperature  distribution,  by 
a  time  of  2  us,  the  dip  is  gone  and  the  temperature  distribution  is  highest 
at  the  center. 

The  current  density  plots.  Fig.  9,  display  some  interesting 
features.  The  conductivity  of  the  air  experiences  a  sharp  cutoff  as  the 
temperature  drops  below  about  2  eV.  Therefore  the  current  is  confined  to 
the  region  with  a  temperature  above  2  eV.  The  current  "spike"  that  occurs 
near  the  outside  of  the  hot  channel  is  caused  by  the  small  skin  depth  of  the 
hot  channel  gas.  As  the  time  derivative  of  the  current  decreases,  the  skin 
depth  of  the  current  increases  and  the  outer  layer  of  increased  current 
thickens,  as  is  seen  in  the  plots. 

The  density  and  pressure  plots.  Figs.  10  and  11,  show  a  shock  wave 
propagating  outward  at  the  edge  of  the  channel.  At  longer  times  than  were 
used  in  the  calculations,  the  shock  wave  separates  from  the  channel  of  hot 
gas  and  travels  outward  through  the  cool  air. 
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The  calculated  optical  emission  of  the  channel  was  obtained  by  con¬ 
volving  the  detector  response  curve  of  the  Lite  Mike  detector  with  the  cal¬ 
culated  frequency  spectrum  of  the  radiated  light.  It  would  be  prohibitively 
expensive  to  use  enough  frequency  groups  to  calculate  the  frequency  spectrum 
of  the  radiated  light  with  adequate  resolution  for  convolution  with  the 
detector  sensitivity  curve.  Therefore,  the  radiated  light  was  assumed  to 
have  a  Planck  distribution  at  the  temperature  of  the  radiating  cell,  within 
the  group  that  contained  the  detector  response  curve.  This  group  covered 
the  1.0  to  3.0  eV  frequency  interval.  The  calculated  emission  and  the 
measured  emission  are  shown  in  Figure  12.  The  calculated  emission  rises 
to  a  peak  of  11  MW/m,  then  falls  to  about  4.5  MW/m  while- the  current  remains 
almost  constant.  The  measured  optical  emission,  however,  rises  to  about 
3.7  MW/m  and  then  remains  constant  for  the  1.5  us  shown  on  the  plot.  This 
discrepancy  probably  results  from  the  initial  conditions  selected  for  the 
calculation.  The  CHARTB  computer  code  does  not  contain  an  electrical  break¬ 
down  model,  but  begins  the  calculation  with  an  initial  small  hot  channel. 

The  size  and  temperature  of  this  initial  channel  affects  the  radiated  output 
for  times  of  a  few  microseconds.  An  interesting  area  for  future  research 
would  be  the  application  of  gas  breakdown  theory  to  the  determination  of 
appropriate  initial  conditions  for  MHD  programs,  such  as  CHARTB. 

The  light  output  from  a  second  ligntniny  stroke  was  also  calculated. 

1  he  measured  current  is  shown  in  Figure  13  and  the  current  waveform  used  in 
the  calculation  "is  shown  in  Figure  14.  The  calculated  optical  emission  con¬ 
volved  with  the  detector  sensitivity  and  the  measured  optical  emission  are 
shown  in  Figure  15.  The  calculated  output  at  a  time  of  1  us  is  about  5  MW/m, 
while  the  measured  value  is  5  MW/m.  This  is  very  good  agreement.  The  spike 
at  the  origin  of  Figure  16  is  probably  dependent  on  the  initial  channel 
temperature  distribution  used  in  the  calculations. 


TIME  (s) 


Figure  14.  Current  waveform  used  to  simulate 


RADIANCE 


4.  RELATION  OF  CHANNEL  LUMINOSITY  TO  IRRADIANCE  AT  LITE-MIKE 


4.1  INTRODUCTION 

Optical  observations  are  used  to  determine  the  illuainance 
(in  lm/m2),  or  equivalently  the  irradiance  (in  W/m2).  The  subject  of  this 
section  is  how  to  obtain  from  this  information  the  channel  luminosity  per 
unit  length  (in  W/m).  This  is  done  as  two  separate  problems:  (1)  relating 
the  channel's  optical  properties  to  the  luminance  and  then  the  brightness 
at  the  film  (or  detector),  and  (2)  relating  the  luminance  to  the  luminosity 
or  radiation  per  unit  length.  We  will  show  that  the  relationships  between 
channel  luminosity  and  observed  brightness  for  optically  thick  and  thin 
channels  differ  only  slightly.  We  will  develop  the  theory  by  treating  pro¬ 
gressively  more  complex  radiators:  a  plane,  Lambert's  law  radiator,  a 
plane,  non-Lamberti an  surface,  a  cylindrical  Lambert’s  law  radiator,  and 
lastly  an  optically  thin  cylinder. 

4.2  THEORY 


4.2.1  Simple  Plane  Radiator 

3ecause  we  nave  found  inconsistencies  between  defi  melons  of  Ijmi  nance 
between  references  (Refs.  7-g)  and  oecause  some  derivations  do  not  list 
underlying  assumptions  (Lambert's  law),  we  feel  it  is  useful  to  define 
terms  and  to  treat  first  the  simple  case  of  the  image  bngntness  of  a  plane 
surface  of  infinitesimal  area,  perpendi  cul  ar  to  the  1  i  rre-of-si  ght .  We  con¬ 
sider  first  a  surface  satisfying  Lambert's  law,  such  as  a  black  body  or 
optically  thick  radiator.  Then  we  will  consider  the  effect  of  relaxing 
this  assumption.  Finally,  we  will  turn  to  the  geometry  of  interest,  that 
of  an  infinite  cylinder. 


We  will  use  SI  jnits  for  radiometric  quantities,  'kite  mat  rao-ant 
iterance  was  romerly  called  radiance  and  radiant  incidence  was  called 
irradiance.  To  convert  to  pnotometric  quantities,  one  would  use  lumens 
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instead  of  watts  in  all  quantities;  radiant  flux  would  become  luminous 
flux,  radiant  incidence  (irradiance)  becomes  illuminance,  radiant  sterance 
(radiance)  becomes  luminance,  etc.  Knowing  the  spectral  distribution  of 
the  radiation,  one  can  convert  from  one  system  to  the  other  using  the  stan-‘ 
dard  luminous  efficiency  curve  (e.g..  Ref.  9  graph,  p.  39,  table  p.  40). 
Table  1  describes  our  notation.  As  we  do  not  consider  any  electromagnetic 
variables,  no  confusion  should  arise  about  the  variable  used.  The  vari¬ 
ables  are  listed  in  order  of  appearance  in  the  text. 

Larmore  (Ref.  8)  finds  the  radiant  incidence  at  the  lens  (retaining 
his  notation  but  changing  from  his  photometric  units  to  radiometric  units) 

I  =  EA/(ir  R2),  where  E  is  the  radiant  exitance  (formerly  emittance  or  raai- 
ant  flux  density)  at  the  source  (in  W/m2),  A  is  the  surface  area  (m2),  and 
R  (m)  is  the  distance  between  the  emitter  and  the  lens  (Fig.  16).  We  may 
relate  the  radiant  incidence  E  at  the  source  and  the  radiant  sterance  L  in 
(W/m2/sr),  in  the  notation  of  (Ref.  9),  by  integrating  over  the  angle  3 
Detween  the  surface  normal  and  the  direction  of  emission  (line  of  sight). 

If  Lambert's  law  holds,  the  radiant  intensity  (W/sr)  is  proportional  to 
cos(0)  and  the  radiant  sterance  L  is  independent  of  9.  Then 

*/2 

E  =  2  it  /  L(  9=0)  cos  (  9 )  sin(9)  d0. 

0 

=  TT  L  (11) 

Page  167  of  Ref.  9  gives  as  the  radiant  flux  F  at  the  lens  F  *  it  LA 
[sin(q)]2,  where  q  is  the  half-angle  subtended  by  the  camera  entrance  pupil 
at  the  emitter.  We  should  have  I  =  F/[n(D/2) 2],  where  D  is  the  lens  diame¬ 
ter.  For  D  «  R,  we  have  sin(q)  =  Q/2R.  We  find  I  =  LA/R2  which  indeed 
checks  with  the  results  of  Ref.  2. 


Variable 
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Radiant  incidence  (irradiance) 
(at  entrance  pupil) 


Radiant  intensity 


Radiant  exitance 


Radiant  flux 


Radiant  incidence  (image) 


Volume  emissivity 


t,e,h  are  dimensionless  factors 
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Figure  16.  Geometry  for  illumination  of  detector  by  planar  radiator. 

We  now  relate  the  flux  or  flux  density  (exitance)  at  the  lens  to  the 
image  brightness  at  the  film.  This  is  i  *  F/a  where  a  is  the  area  of  the 
image.  We  have  i  =  I ttD 2/4a  =  EAD2/4aR2.  As  the  magnification  m  is  related 
to  the  areas  A  and  a  and  the  optical  system's  focal  lenyth  f  by  m2  =  a/ A  = 
(f/R)2  we  have  i  =  ED2/4f2,  as  given  by  Larmore  (Ref.  8).  Note  that  the 
results  are  independent  of  the  distance  between  the  source  and  the  film. 
Note  also  that  we  have  neglected  transmission  losses,  both  within  the  opti¬ 
cal  system  as  well  as  between  the  source  and  the  camera.  These  losses  will 
be  accounted  for  by  a  simple  multiplier. 

For  a  plane  radiator,  we  are  concerned  with  a  radiative  loss  rate 
E  in  (W/m2)  but  not  a  loss  per  unit  length.  In  this  case  we  have  simply 
'•elated  i  and  E  and  no  more  work  need  be  done.  We  might  write  E  =  ii,t  (f/ 
wnere  t  is  the  net  transmission  ractor  for  the  atmosphere  and  f;D  is  the  f- 
of  tne  optical  system.  This  would  solve  our  task  of  relatina  E  to  i . 


4.3  NON-LAMBERTIAN  RADIATOR:  OPTICAILY-THIN  PLANE  SLAB 


If  the  plane  emitter  is  not  a  blackbody  source,  but,  e.g.,  an  opti¬ 
cally  thin  layer  of  emitting  gas,  we  would  find  in  place  of  Lambert's  law 
(L  independent  of  9),  L  *  sec  0  due  to  path  length.  The  radiant  inten¬ 
sity  would  be  independent  of  8,  as  the  factors  of  sec  9  and  cos  9  would 
cancel.  Then  we  have  L(9=0)  =  E/2ir,  and  ultimately  E  =  8i/t(f/D)  2.  This 
is  a  factor  of  two  different  from  our  previous  results.  Finite  optical  - 
depth  sources  will  be  intermediate  between  the  two  cases  presented  here,  a 
factor  of  the  form  1  +  exp(-d)  where  d  is  the  optical  thickness  of  the 

slab,  may  be  used  as  an  approximate  interpolation  factor.  Note  that, 

because  of  the  small  angle  approximation,  the  departure  from  Lambert's  law 
does  not  materially  affect  the  relationship  between  i  and  I  or  between  I 
and  L,  only  the  relationship  between  E  and  L. 

4.4  INFINITE  CYLINDER,  LAMBERT'S  LAW  RADIATOR 

We  now  proceed  to  develop  similar  results  for  an  infinite  cylinder. 
Applying  the  formalism  of  Ref.  9,  we  have: 

l  it/ 2  2ir  q 

F  =  2 r  /  dz  /  de  /  d$  /  [cos  p  L(i>)]  sin  (a)  da 

0  U  0  0 

where  the  angles  are  defined  in  Figure  17.  The  small  angle  approximations 
simplify  things  considerably,  as  the  cosine  of  any  small  angle  may  oe  taken 
as  1,  to  first  order.  Then  p  =  r/R  sin  8  -  sin  a  =  (1  +  r/R)S  -  a  =  9  to 
the  required  accuracy.  Consider  a  section  of  the  channel  of  length  1  for 
concreteness.  We  then  find  F  =  wL  2 ri  sin2(q),  i.e.,  the  the  same  as  for  a 
plane  surface  emitter  with  the  same  projected  area  (area  perpenddcul ar  to 
the  line  of  sight),  as  we  might  have  expected.  Setting  the  projected  chan¬ 
nel  area  A  =  2 rl,  we  nave  again  I  =  LA/R2  and  i  =  L  it/ 4  ( □/ f )  2 ,  but  now 
instead  of  L  =  E/t  we  nave  to  relate  L  to  a  loss  rate  per  unit  c.nannel 
iengtn.  For  an  optically  tnicx  channel  we  nave  A  {'A/m)  =  2u  rS  wnere  the 
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constant  (~1012  for  T  in  eV),  and  L  =  S/w  since  the  black  body  emitter 
obeys  Lambert's  law.  Thus,  W  =  2n2rL,  and,  inserting  the  transmission  fac¬ 
tor  t, 

U  =  h2R2I/  ( it)  =  8irri/t  (f/D)2.  (12) 

This  is  the  required  relation  for  channel  emission  W  from  the  data  i  and 
r.  The  latter  may  be  found  from  the  image  size  s,  knowing  the  camera  sys¬ 
tem  magnification,  sr  =  mr  and  m  =  f/R. 

4.5  OPTICALLY  THIN  CYLINDRICAL  CHANNEL 

As  aDove,  we  must  find  a  new  relation  between  W  and  L(0),  as  well  as  a 
relation  between  i  and  L(0).  In  a  plane  perpendicular  to  the  channel  axis, 
we  have  an  L(0)  proportional  to  cos  9-;  in  a  plane  containing  the  channel 
axis,  we  have  L(9)  proportional  to  sec  9.  Note  that  the  film  image,  if  it 
were  to  possess  sufficient  resolution,  contains  information  about  the  opti¬ 
cal  depth  of  the  channel;  if  the  channel  limb-darkening  corresponds  to  a 
cos( 9 )  law,  we  may  take  the  channel  as  thin.  If  the  brightness  seems  con¬ 
stant  across  the  channel  cross-section,  we  may  treat  it  as  optically 
thick.  Unfortunately,  in  typical  cases  films  have  one  to  three  grains 
across  the  channel.  Consequently,  theory  must  be  used  to  estimate  the  , 
optical  depths  in  tne  spectral  re'jime  of  interest  -film  response).  "hat 
will  not  oe  treated  oere. 

Assuming  the  cnannel  is  optically  thin,  we  acquire  ^per  unit  lenytn), 
to  a  jood  approximation,  an  addition  factor  of  cos( 9)  in  the  integral 
relating  F  to  L { 0 )  .  The  integral  over  9  is  now  n/4,  not  1  as  in  the  opti¬ 
cally  thick  case.  This  is  a  decrease  of  only  21%.  We  have  in  general 
I  =  eL(Q)A/R4  where  e  =  1  in  the  truck  limit  and  tt/ 4  in  the  tmn  limit. 

For  tne  tmn  case,  tne  loss  rate  oer  unit  length  W  (W/m)  =  e  n  r‘  wnere  £ 
is  tne  vol.ne  emissiv'ty  m  watts  oer  meter3.  We  know  that  L;-j)  =  _'r  c.ot, 
oecause  tie  outi  til  Path  '  engt.n  is  2r  'or  v  -  ->  -  ij,  and  each  volume  ele¬ 
ment  along  vi  ,jt  jd  tn  ra<j  i  ales  ■  so  L  rup  i  .a  i  1  y  over  4u  steraoians. 
Consequently,  a  volume  element  dV  contnoutes  £dV/4ir  Watts/ steradi  an  to  the 
radiant  intensity  in  any  direction.  Note  tnat  as  the  radiant  intensity 


contribution  of  any  element  is  independent  of  the  distance  from  that  ele¬ 
ment,  we  may  simply  integrate  along  the  line  of  siyht  to  obtain  the  radiant 
sterance  in  any  direction.  Thus,  in  the  optically  thin  case, 

W  =  L(Q)  r2tr2,  as  in  the  thick  case.  Then 

W  *  w2IR2/(ei)  .  (13) 

Finally  i  =  ne/4  L(0)  (D/f)2  yielding  the  answer  W(Watts/m)  =  i(Watts/n2) 
x(f/D)2  8irR/(et).  Given  I  or  i ,  there  is  at  most  a  21%  error  in  estimating 
W,  given  uncertainties  in  the  channel  optical  depth. 


The  power  incident  at  a  detector  of  area  A  is  P  =  AI.  Thus,  in 
general , 


where 


e  =  1  thicK  channel 
4/ n  thin  channel 


or 


q  =  1  thin  channel 

4 / *  thick  channel . 


This  result  (q  =  1)  is  used  above  in  the  analysis  of  Lite  Mike  data. 

In  many  applications  the  result  of  a  microdensitometer  scan  across  the 
cnannel  will  oe  the  observable:  B  (W/m)  =  j  i  dx.  For  an  optically  tnicx 
channel  there  is  no  limb  darkening  and  B  =  iw  where  w  is  the  image  width  of 
the  cnannel.  For  a  thin  channel,  while  the  image  brightness  will  vary  as 
cos ( 0 ) ,  if  we  define  w  similarly  we  will  have  the  same  relationship  between 
B  and  i.  This  is  because  luminance  is  conserved  (except,  of  course,  for 
the  losses  accounted  for  by  the  transmission  factor  t)  and  the  relationship 
oetween  image  flux  density  and  entrance  pupil  flux  density  and  entrance 
pupil  flux  is  unchanged  (Ref.  9).  One  may  also  argue  rrom  a  conservat i on 
of  energy  viewpoint.  Thus,  we  find 
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„  f2  8nBh  R 
w  =  TT  ~t - w 

where  h  =  1  if  the  channel  is  thick  and  4/ir  if  thin.  But  2R/w  =  p/f  =  1/m 
(see  Section  4.4),  so  that 

w  s  4iih  f£_  g  .  (14) 

t  D2 

Note  W  may  be  inferred  from  B  without  knowledge  of  the  channel  radius,  with 
channel  optical  depth  introducing  only  a  small  uncertainty  if  unknown. 

The  optical  depth  of  the  channel  must,  however,  be  Huiomi  to  find  h, 
and  the  spectral  character! sties  of  the  radiation  must  be  known  to  relate 
the  radiant  quantities  (i.e.,  those  specified  in  Watts)  to  the  film 
response. 

4.6  SIMPLER  DERIVATION  BASED  UPON  SYMMETRY 


.re 


We  assume  that  a  section  of  a  lightning  channel  has  been  photographed 
and  that  the  section  of  channel  posesses  axial  symmetry  with  the  axis  per¬ 
pendicular  to  the  line  from  the  section  of  channel  to  the  camera  lens. 


We  snail  obtain  the  output  power  passing  through  the  surface  of  a  sec¬ 
tion  of  length  I  of  the  channel  by  integrating  the  corresponding  irraoiance 
over  the  surface  of  a  sphere  of  radius  R,  that  is  centered  at  the  section 
of  channel.  The  camera  aperture  diameter,  D,  is  small  compared  to  R. 
Therefore  I,  the  irradiance  (total  flux  per  unit  area),  is  uniform  across 
the  camera  lens,  and  corresponds  to  the  irradiance  at  90°  to  the  channel 
axis.  Consider  a  detector  measuring  the  irradiance  at  a  point  on  the  sur¬ 
face  of  the  sphere.  Let  9  and  <j>  be  the  coordinates  of  the  point  on  the 
spnere.  The  angle  the  line  from  the  center  of  the  section  of  channel  to 
the  detector  makes  with  the  channe1  axis  is  9,  and  the  angle  about  the 
channel  axis  is  f>.  The  irradiance,  I,  is  independent  of  ?  but  will,  in 
general,  depend  on  9.  [  is  the  radiant  incidence  or  irradiance;  it  is  the 

flux  per  unit  area  at  a  point  on  the  sphere  that  come  from  all  locations 
within  the  lightning  channel.  We  can  limit  our  discussion  to  Il(  t^e 
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irradiance  passing  through  the  surface  of  a  short  length  (z  «  R)  of  the 
channel,  by  considering  only  the  photographic  image  of  that  portion  of  the 
channel . 


For  visible  wavelengths,  the  channel  is  optically  thin.  The  irradi¬ 
ance  from  the  length  of  channel  is  proportional  to  the  solid  angle  sub¬ 
tended  by  the  length  of  channel  as  seen  from  the  viewing  location, 

I£  «  sin  0. 

For  an  optically  thin  channel,  the  irradiance  is  proportional  to  the  aver¬ 
age  path  length,  within  the  section  of  channel  of  rays  arriving  at  the 
viewing  location, 

I  <x  ,  1 

sin  S  * 


Therefore,  for  an  optically  thin  channel,  1^(9)  is  independent  of  9, 
Iz(0)  =  I4(90°) 


and  W,  the  power  output  of  the  channel  in  watts  per  meter  may  be  obtained 
by  integrating  I4  over  the  surface  of  a  sphere  of  radius  R. 


3  2  tt/2  2tt 

W  =  2  f-  /  d9  /  do  1,(90°)  sin  9 
0  0 


W 


4ttR2 

Z 


I£(90°) 


ir/2 

r 

i 

J 

0 


d9  sin  9  = 


4ttR2 

z 


Equations  16  and  13  agree. 


(15) 


(16) 


The  irradiance  at  the  camera  lens  is  obtained  by  integrating  If,  the 
irradiance  at  the  film  over  the  image  of  the  channel  section. 


where  M  is  the  image  magnification,  y  is  the  direction  along  the  channel,  x 
is  across  the  channel,  and  A  is  the  area  of  the  camera  aperture.  The 
assumption  that  If  is  independent  of  y  leads  to  the  simplification 


.  MA 

h  s  T 


/  U  dx  . 


The  integral  immediately  above  may  be  evaluated  even  if  the  channel  is  not 
wel 1 -resol ved  due  to  the  graininess  of  the  photograph.  The  equation  for 
the  power  output  per  unit  channel  length  is 


4*R2M 


/  If  d'x 


The  integration  limit,  a,  must  be  taken  large  enough  to  include  the  width 
of  the  channel  on  the  photograph. 

Introducing  N  =  the  camera  f-number,  F  =  the  camera  focal  length,  and 
D  =  the  camera  aperture  diameter,  the  camera  aperture  area  is 

'  A  3  j  s  \  ^N2F2 

and  the  magni fication  is 
M  «  F/R  . 

The  channel  radiated  power  per  unit  length  is 


rF  -a 


If dx . 


Equation  20  is  our  main  result. 

It  may  oe  of  interest  to  obtain  W  for  an  optically  thick  channel.  The 
the  irradiance  at  an  observation  point  originating  at  the  surface  of  the 
segment  of  channel  is 

I a( 9)  =  Iz{y0°)  sin  d  f  21 ) 
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because  the  intensity  of  a  ray  is  independent  of  the  distance  it  has  passed 
through  the  hot  gas.  The  power  output  of  the  channel  is 

R2  V2  2 

W  «  2  7-  /  d8  /  d*  I .(90°) sin 28 

1  0  0  4 
U  *  1,(90°)  d8  \  Z  sin2  8 

i.  i.  o 

2  2 

W  =  ^ -  I£(90°)  .  (22) 

Equations  26  and  12  agree.  Therefore,  for  an  optically  thick  channel,  we 
use  Equations  18  and  22  to  obtain 

W  *  /3  If  dx  (23) 

or 

W  =  M  /3  i  dX  (24) 

N2F  -a 

where  W  is  in  watts  per  meter.  Equation  24  gives  the  power  output  for  an 
optically  thick  channel,  and  is  a  factor  of  (ir/4)  times  the  power  radiated 
by  an  optically  thin  channel  (Eq.  20).  In  applying  Equations  20  or  24,  it 
is  not  necessary  to  determine  the  width  of  the  channel,  nor  is  it  necessary 
to  resolve  the  brightness  at  a  specific  point  on  the  channel. 

4.7  CONCLUSION 

The  principal  conclusion  is  that  optical  observations  of  lightning  can 
yield  accurate  energy  loss  rates  without  knowledge  of  channel  radius  or 
optical  depth.  The  formulae  above  are  ready  for  application.  Only  the 
optical  depth  in  the  spectral  region  of  interest  is  needed,  as  well  as 
estimates  of  the  channel  radius  (uncertain  due  to  film  resolution  by  a  fac¬ 
tor  of  roughly  two  to  three).  Note  that  Equation  20  or  24  may  be  applied 
to  a  case  in  which  the  channel  radius  is  not  resolved,  without  the 


necessity  of  estimating  the  channel  radius.  The  parameter  E  may  be  inter 
polatedas  a  function  of  optical  depth  for  the  cylindrical  case  as  in  the 
slab  case;  optical  depths  of  a  few  are  indistinguishable  fro*  infinite 
optical  depth,  depths  smaller  than  one  are  qualitatively  similar  to  the 
thin  limit.  Even  a  poor  estimate  of  the  channel  optical  depth  will  yield 
an  accurate  loss  rate.  The  atmospheric  transmission  is  another 
uncertainty. 


5.  ANALYSIS  AND  INTERPRETATION  OF  B-DOT  SIGNALS 


5.1  RESULTS  OF  COMPUTER  MODELING  OF  FIELDS  FROM  THE  STROKE  CURRENT 

The  Lite  Mike  experiments  provided  data  from  B-Dot  and  D-Dot  sensors 
located  at  Kiva  1,  in  addition  to  the  Lite  Mike  and  current  sensor  data. 

It  was  felt  that  it  would  be  of  interest  to  analyze  these  data  records. 
While  a  great  deal  of  data  from  these  sensors  have  been  reported  in  Lightn¬ 
ing  Phenomenology  Note  3  and  Lightning  Phenomenology  Note  11,  a  detailed 
analysis  of  these  data  does  not  seem  to  have  been  performed  for  any  cases 
yet.  We  hope  that  the  results  presented  here  stimulate  the  application  of 
sophisticated  models  to  the  treasury  of  data  from  the  Kivas.  One  novel 
discovery,  that  of  a  higher  initial  return  stroke  velocity  than  previously 
observed,  has  been  made.  It  is  necessary  to  confirm  this  result  with  addi¬ 
tional  observed  cases.  The  confirming  data  may  already  exist,  and  merely 
await  detailed  analysis. 

The  strike  of  record  85220  will  be  analyzed.  It  appears  to  be  rather 
typical  of  strikes.  There  are  two  strikes  for  85222;  these  posses  anoma¬ 
lously  low  current  amplitudes  of  about  2000  A  maximum.  Furthermore,  the 
8-Dot  and  D-Dot  records' have  not  been  made  available  for  analysis. 

The  measured  current  of  strike  85220  is  shown  in  Fig.  6.  The 
polarity  is  positive  vertical  current  (i.e.,  the  lowering  of  negative 
charge  or  the  rise  of  positive  charge),  which  is  the  most  common  form  for 
rocket-tri ggered  lightning.  This  current  was  modeled  for  computational 
purposes  as  a  double  exponential  waveform,  shown  in  Fig.  18.  The  fall  time 
is  not  important  as  the  results  are  most  sensitive  to  the  rise  time. 

Videotape  observations  show  that  the  stroke  was  approximately  verti¬ 
cal  within  the  vicinity  of  the  stinger,  witn  an  inclination  of  approxi¬ 
mately  15  degrees  to  the  vertical,  the  channel  moving  toward  the  soutneast 
with  increasing  height.  We  will  consider  a  vertical  channel,  and  treat  the 
inclination  of  the  channel  through  an  approximation.  The  vertical 
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Figure  18.  Measured  and  computed  current  irr  chairomel 
versus  time. 


component  of  the  current  will  be  used  to  compute  the  vertical!  electric  dis¬ 
placement  and  the  azimuthal  magnetic  induction  at  the  location  of  the  sen¬ 
sors,  and  the  distance  will  not  be  measured  from  the  vertical  line  but  from 
the  displaced  channel  location,  taken  to  be  in  the  EU  plane. 


The  uasic  method  is  that  of  liman,  McLain,  and  Krider  (Ref.  ID),  who 
consider  an  ideal  (lossless)  transmission  line  with  the  current  waveform 
propagating  upward  without  modification.  The  ground  is  assumed  to  be  per¬ 
fectly  conducting,  producing  ideal  image  currents  in  the  ground  plane. 
Because  of  the  conducting  wire  grid  in  the  area  of  the  sttnger,  this  should 
be  an  excellent  approximation.  The  ground  is  not  perfectly  flat  or  uni¬ 
form,  with  the  sensors  at  Kiva  1  approximately  4  m  below  Kiva  Z.  The  sen¬ 
sors  will  oe  assumed  to  oe  4U  m  from  the  channel  base.  The  stinger  neignt 
is  20  m.  To  take  full  account  of  the  physics  of  stroke  initiation,  one 
would  have  to  consider  not  one  wave  travelling  upward  from  the  ground  (ana 


its  image),  but  one  wave  starting  upward  from  the  initiation  point,  and 
another  starting  downward  from  this  point.  These  waves  would  have  images. 

In  addition,  the  downward  wave  would  be  reflected  from  the  ground,  giving 
rise  to  an  upward  wave.  To  a  good  approximation,  these  added  complications 
may  be  ignored  after  the  transit  of  the  wave  back  to  the  initiation  point, 
i.e.,  after  a  time  2H/c,  where  H  is  somewhat  greater  than  20  m  and  c  is  the 
speed  of  light.  This  time  is  roughly  130  ns.  We  will  make  this  approxima¬ 
tion,  which  is  justified  a  posteriori  by  the  agreement  of  the  observations 
and  computational  results.  The  wave  is  assumed  to  propagate  at  a  velocity 
V. 

The  azimuthal  magnetic  induction  B  is  the  sum  of  two  terms,  an  induc¬ 
tion  term  and  a  radiation  term,  the  former  having  an  inverse  square  depen¬ 
dence  upon  distance  from  the  radiating  element,  and  the  latter  an  inverse 
dependence  upon  distance.  This  results  in  the  radiation  term  dominating  at 
very  large  distances.  At  the  very  close  distances  of  interest  here,  the 
terms  are  of  comparable  magnitude.  At  very  early  times,  the  terms  are  both 
positive  with  the  radiation  term  dominating  due  to  the  rapidity  of  the 
field  changes.  After  400  ns  in  the  case  computed  and  displayed  in  Figures 
19  through  22,  the  induction  term  is  larger  than  the  radiation  term.  At  late 
times  (after  1.7  us  for  the  V  =  0.9  c  wave  velocity  case  considered  here) 
the  radiation  term  becomes  negative,  but  it  typically  remains  an  order  of 
magnitude  smaller  than  the  induction  term.  The  full  expression  for  3  is 
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(25) 


wnere  e  is  the  angle  between  the  channel  and  the  radius  between  the  radi¬ 
ating  element  and  the  observation  point  on  the  ground  plane.  Thus,  each 
term  is  positive  and  decreases  rapidly  as  we  go  up  the  channel,  ooth 
aecause  of  tne  ‘actors  of  R  and  the  factor  of  sin!  el. 
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Figure  21.  Observed  and  computed  magnetic  induction  B 
versus  time  for  wave  velocity  0.9  c. 


The  electric  field  is  more  complicated,  having  a  static  term,  an 
induction  term  and  a  radiation  term: 
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Figure  22.  Observed  and  computed  derivative  of  magnetic 
induction,  B. 


The  static  term  is  due  to  the  effective  dipole  moment  formed  by  the 
moving  current  pulse  and  it  image,  which  results  in  a  1/ R 3  dependence. 
This  term  does  not  vanish  at  large  times  due  to  the  finite  charge  trans¬ 
ported  oy  the  stroke.  All  the  terms  contributing  to  the  electric  field 


and  electric  displacement  have  much  stronger  dependences  upon  channel 
height,  due  either  to  a  factor  of  sin29  instead  of  sin  in  the  radiation 
term  or  factors  of  2-3  sin29  in  the  other  terms,  and  the  R3  dependence  of 
static  field  term.  Note  that  the  static  and  induction  field  terms  are  of 
indefinite  sign.  Lower  levels  of  the  channel  have  values  of  sine  =  D/R 
close  to  one,  giving  a  negative  contribution.  Higher  portions  of  the  chan¬ 
nel  have  small  values  of  the  sine  term,  giving  a  positive  contribution. 

The  factor  (2-3  s i n 2 8 )  does  not  vanish  with  height  as  do  the  factors  of  sin 
in  the  terms  for  B,  so  the  induction  term  might  be  more  significant  than 
the  correspondi ng  term  in  the  equation  for  B.  The  crossover  occurs  at  a 
height  of  1.414  0,  where  D  is  the  distance  of  the  observer  from  the  channel 
base.  In  our  ca-se  this  is  a  mere  60  m  or  so.  It.  is  seen  from  the  downward 
nature  of  D  and  E  for  upward  current  that  the  portions  of  the  channel  below 
this  height  dominate  in  their  contributions  in  the  measured  signal. 

However,  the  static  term  is  found  to  dominate  the  other  terms  except  at 
very  early  times,  at  which  all  play  a  role.  Initially,  all  terms  are  nega¬ 
tive,  with  the  radiation  term  dominant  due  to  the  large  time  derivatives. 
Soon  the  static  term  dominates,  on  the  same  timescale  as  the  induction  term 
in  the  magnetic  induction  dominates  the  radiation  term.  Both  the  induction 
and  radiation  terms  change  sign,  the  former  at  1  ns  and  the  latter  at 
1.4  us  in  the  case  of  Figs.  19  through  22.  The  static  term  is  about  one 
order  of  magnitude  larger  than  the  induction  field  which  in  turn  is  one 
order  of  magnitude  larger  than  the  radiation  term.  Obviously,  these 
rel ati onshi ps  will  change  if  the  observer  is  located  farther  from  the  the 
channel  or  if  the  current  pulse  differs. 

For  computational  purposes,  the  radiation  terms  were  integrated  by 
replacing  the  time  derivative  of  the  current  with  the  derivative  of  current 
with  respect  to  height.  Thus,  the  second  integral  in  the  expression  for  3 
becomes 


J  J i  ■  si nu 
2mC  '  A 


di 


■15 


....  /■/%■ . 


Y«\V.  .v J-f-r.  ■',  •  . 

r  •  -  *  **  *  "_»  r  t 


.  v.v- 

&  —  *  ii.A  it. £  L t . 


which  may  be  integrated  by  parts  to  yield 
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The  radiation  term  in  the  equation  for  the  electric  field  is  similar  except 
that  sin^e  replaces  sine.  With  the  substitution  sine  *  D/R  and 

- * - y— 

R  =  /(l)"  +  z  )  we  may  calculate  the  derivative  in  the  integral,  which  is 
-2z/(D2  *  z2)  for  the  expression  above. 
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The  static  field  term  is  calculated  by  defining  a  variable  q(z)  as 
the  time  integral  of  the  current  at  any  height.  As  the  observer  is  sta¬ 
tionary,  this  suffices  because  the  retardation  to  any  point  on  the  channel 
does  not  change. 

In  the  transmission-line  model  the  principal  free  parameter  (once  the 
waveform  of  the  current  pulse,  which  is  completely  free,  is  specified),  is 
the  wave  velocity.  Good  agreement  was  found  for  the  model  used  with  a  wave 
velocity  of  V  =  G.9  c,  where  c  is  the  speed  of  light.  Ue  present  the 
results  for  this  assumed  wave  velocity  as  follows:  Fig.  19  compares 
observed  and  computed  electric  displacement  D;  Fig.  20  compares  D-Oot ; 

Fig.  21  compares  magnetic  induction  B,  and  Fig.  22  B-Oot.  Note  that  the 
sensor  for  D-Oot  causes  a  sign  change  between  the  computed  and  observed 
results  for  these  quantities. 

Figures  23  through  26  give  results  for  computed  fields  with  an 
assumed  wave  velocity  of  c/3.  It  will  be  seen  that  the  peaks  in  the  deri¬ 
vative  signals  are  delayed  and  reduced  in  magnitude,  resulting  in  poorer 
agreement  with  observations. 

5.2  INTERPRETATION  OF  RESULTS 

The  wave  velocity  which  best  matches  observations  is  higner  than 
typically  reported  for  return  strokes.  Oman  (Ref.  11)  gives  c/3  as  a  typi¬ 
cal  values  and  c/2  as  a  typical  upper  limit  to  reported  velocities. 
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Figure  23.  Computed  results  for  wave  velocity  c/3: 
Electric  displacement  D. 
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Figure  24.  Computed  results  for  wave  velocity  c/3: 
Magnetic  induction  3. 
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Figure  25.  Computed  results  for  wave  velocity  c/3: 

Derivative  of  electric  displacement,  D. 
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Figure  26.  Computed  results  for  wave  velocity.c/3 : 

Derivative  of  magnetic  induction,  3. 


The  explanation  most  likely  is  that  near  the  ground  the  velocities  are 
higher  than  average.  The  wave  is  traveling  in  the  most  recently  formed 
portion  of  the  channel,  which  is  probably  still  highly  conductive  and  wnich 
has  not  dispersed  to  any  degree.  As  the  wave  progresses,  it  encounters 
channel  segments  which  have  cooled  and  undergone  attachment,  recombination 
and  some  hydrodynamic  dispersal  due  to  turbulence,  etc.  The  corona  enve¬ 
lope  which  would  store  the  charge  around  the  channel  is  probably  of  greater 
diameter  as  well  as  more  resistive  due  to  attachment.  This  would  account 
for  the  reduced  wave  speed  (Refs.  12,  13). 

The  high  wave  velocity  is  also  of  great  interest  as  it  results  in  a 
correspondingly  larger  value  of  radiated  fields  for  a  given  current  than 
for  low  wave  velocity.  This  is  of  importance  in  comparing  the  I  of  return 
strokes  inferred  by  different  methods  (e.g.,  remote  and  tower  measure¬ 
ments)  (Ref.  14). 

In  the  present  work,  it  has  been  adequate  to  assume  constant  wave 
speed  and  wave  shape  to  obtain  reasonable  agreement  with  observation. 
Residual  errors  might  be  explained  due  to  the  imperfect  ground  plane 
(neither  flat  nor  perfectly  conducting)  and  the  non-ideal  channel  geometry 
(not  perfectly  vertical).  Because,  as  noted  above,  the  very  lowest  por¬ 
tions  of  the  channel  are  dominating  the  signal,  it  is  quite  possible  tnat 
the  wave  has  moved  beyond  this  region  of  importance  before  it  has  changed 
significantly  in  profile  or  velocity.  As  this  region  of  significance  is 
roughly  the”  lowest  60  m  of  the  channel,  20  m  of  wnich  is  contributed  by  the 
stinger,  this  should  not  be  surprising. 

The  details  of  wave  initiation,  which  would  influence  the  first  10  ns 
or  so  of  the  signal,  do  not  seem  to  be  important  in  contributing  to  the 
observed  fields. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  these  results,  we  would  conclude  the  following: 

(1)  Near  field  values  (B  and  D)  can  be  adequately  explained  with 
simple  models. 

(2)  Wave  velocities  near  the  ground  appear  to  approach  the  speed  of 
light.  As  this  is  significantly  higher  than  reported  wave 
velocities  in  the  literature  for  return  strokes;  this  suggests 
that  the  wave  slows  as  it  rises.  It  is  quite  plausible  that 
the  wave  changes  profile  as  well  as  it  rises.  This  can  be  com¬ 
puted  sel f-consistently  (Ref.  13).  Such  a  computation  is  not 
necessary  for  explaining  the  fields  at  a  sensor  40  m  from  the 
stinger,  as  only  heights  of  order  60  m  or  below  contribute  sig¬ 
nificantly  to  the  fields  here. 

We  would  recommend  for  future  investigations  the  following: 

(1)  Analyses  similar  to  that  performed  here  for  B-Dot  an  D-Dot  data 
as  contained  in  Lightning  Phenomenology  Note  11  for  selected 
strikes.  This  would  give  insignt  into  stroke  initiation 
behavior.  Such  a  series  of  analyses  would  confirm  or  refute 
the  high  initial  wave  velocities  discovered  in  this  work,  and 
would  thereby  contribute  new  facts  to  the  understanding  of  the 
lightning  return  stroke. 

(2)  More  sophisticated  analysis  using  the  transmission-1 ine  model 
of  Ref.  13.  This  would  verify  the  hypothesis  of  the  slowing  of 
the  wave  with  increasing  channel  height.  Most  importantly,  it 
would  enable  us  to  build  a  model  in  which  the  wave  shape  is  not 
arbitrary  out  determined  by  the  physics  of  the  lightning  chan¬ 
nel  . 
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(3)  Locating  other  sensors  at  greater  distances  from  Kiva  2  to 
obtain  information  on  the  wave  character  at  greater  heights. 
This  would  be  useful  only  if  better  information  on  channel 
geometry  were  available  as  well.  This  information  might  in 
part  be  obtained  by  locating  sensors  at  different  a*imuths  from 
the  stinger,  but  probably  would  require  more  videotape  loca¬ 
tions. 
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6.  LIGHTNING  PROTECTION  AND  TRIGGER  INC  fmi'mS 


Lightning  protection  is  usually  accomplished  by  pTrawiiatiing  an  alter¬ 
nate  path  for  the  discharge.  The  alternate  path  mustt  fmrTiute  an  electrode 
that  will  form  a  leader  and  connect  with  any  lightning  Teestaar  before  a 
leader  from  the  protected  object  can  do  so.  The  questtixni)  <xff  srfiether  or  not 
a  leader  from  an  electrode  will  be  formed  and  sustain®*  tty/  a  given  gap  vol¬ 
tage  and  geometry  is  also  of  interest  in  the  analysis  off  Tiigjttning  trig¬ 
gering.  The  structure  of  an  empirical  model  for  ligjttrmhmj  assessment  is 
described  in  Section  6.1.  There  are  unknown  parameters  tJrafc  wist  be  pro¬ 
vided  for  the  model  by  performing  experiments  with  long  ^aarlc  gaps;  or,  it 
may  be  possible  to  deduce  the  parameters  after  conducting;  a  thorough  searcn 
of  the  literature  that  is  available  on  the  breakdown  erf  Tong  spark  gaps. 

The  charging  rate  of  a  grounded  imperfect  conductor  suspended  in  a 
vertical  gap  is  discussed  in  Section  6.2.  This  analysts  ts  related  to  the 
question  of  whether  or  not  the  channel  formed  by  the  detonation  of  a  sus¬ 
pended  detonating  fuse  would  retain  its  conductivity  for  Tong  enough  for 
charge  to  accumulate  at  the  upper  end  so  that  an  upward]  reader  could  be 
launched. 

6.1  #  EMPIRICAL  MODEL  FOR  LIGHTNING  PROTECTION 

A  model  for  lightning  protection  assessment  is  concerned  with  posi¬ 
tive  streamers  that  are  launched  by  objects  near  the  groumij.  In  experi¬ 
ments  conducted  with  long  sparks,  a  given  voltage  waveform'  will  have  a  pro¬ 
bability  of  causing  a  breakdown  of  the  gap.  There  is  a Ts®  a  statistical 
distribution  of  the  time  to  leader  inception  and  the  time  to  breakdown. 

For  protecting  high-voltage  transmission  lines  from  breakdown ,  it  is  only 
important  to  treat  the  probability  of  breakdown  of  a  give??  gap.  Semi- 
empirical  mocels  are  well  developed  for  tnis  case.  For  the  case  of  liynt- 
ning  protection,  however,  the  voltage  at  the  lightning  rod  and  the  pro¬ 
tected  conductor  continues  to  increase  as  the  cloud  to  earth  leader 
approaches  until  an  upward  leader  from  one  of  them  reach®*  it  and  causes  a 


conducting  path  to  ground.  Therefore,  the  probability  distribution  of  the 
relative  times  at  which  leader  inception  occurs  at  the  two  conductors  is  a 
crucial  consideration  for  the  lightning-protection  case,  and  the  existing 
models  can  not  be  used. 


Reference  15  includes  a  brief  discussion  of  the  time  to  breakdown  of 
air  gaps.  We  assume  a  configuration  consisting  of  a  lightning  rod  of 
height  hr,  a  protected  object  of  height  hp,  and  a  negative  leader 
approaching  with  speed  v  from  above  at  height  h4  at  time  t  =  0.  Usually, 
we  expect  that  hA  >>  hr  >  hp.  A  positive  leader  will  be  launched  by 
the  object  or  by  the  lightning  rod,  or  by  both.  At  least  one  of  the  posi¬ 
tive  leaders  will  join  the  negative  leader,  determining  where  the  lightning 
will  strike.  To  illustrate  the  general  method  proposed,  it  is  assumed  that 
the  voltage  between  the  positive  leader  and  the  lightning  rod  is,  to  first 
order, 


ur(t) 


Ui 

dL  -  vt 


where  d^  *  h4  -  hr.  Similarly,  the  voltage  between  the  positive 
leader  and  the  object  is 


up(t)  ■  1 

where  d2  =  hj,  -  hp.  The  time  to  breakdown  of  air  gaps  can  be  fitted  to 
a  Weibull  distribution  for  the  probability  density  (Ref.  15). 

p(y)  =  (o/t)  ya_1  exp(-ya)  ( 


where  y  =  (t-sl/r,  typically,  a  is  given  the  value  a  =  1.5,  and  e  repre¬ 
sents  a  minimum  time  for  breakdown,  e  and  t  are  functions  of  the  breakdown 
overvoltage  factor 


The  probability  of  breakdown  between  the  positive  leader  and  the  liyhtniny 
rod,  during  an  elapsed  time  t  is 

P  l(t)  =  /  n[t'-ei(t')]  [  <*/ 1 !  ( t ' )  ]  yi*”1  exp(-y1a)dt‘  (35) 

0 

where  n  is  the  unit  step  function 


1  for  z  >  0 

n(z)  *  (36) 

0  for  z  <  0 

and  Tj,  e\,  and  yt  are  obtained  by  substituting  Equation  36  into  Equations 
31,  32,  33  and  34.  A  similar  equation  may  be  obtained  for  ?2(t),  the 
probability  that  breakdown  occurs  between  the  positive  leader  and  the  pro¬ 
tected  object  before  time  t. 


This  is  only  the  outline  of  a  possible  semi -empi ri cal  model  for 
lightning  protection.  Verification  of  the  forms  of  the  equations  and  the 
evaluation  of  e0,  e^,  ts,  3,  and  p  for  various  electrode  geometries 
would  require  a  detailed  survey  of  the  available  data  on  breakdown  condi¬ 
tions  for  long  sparks.  More  than  one  parameter,  e.g.,  pulse  rise  time  and 
peak  voltage,  are  frequently  varied  from  measurement  to  measurement,  so  the 
dependence  on  a  single  variable  is  frequently  difficult  to  extract.  The 
type  of  measurement  that  would  be  the  most  useful  for  lighting  protection 
problems  is  not  directly  applicable  to  the  breakdown  of  electrical  power 
transmission  networks,  and  is  therefore  rarely  made.  This  is  the 


measurement  of  the  time  to  breakdown  as  the  voltage  is  increased  at  an 
approximately  constant  rate.  The  time  to  breakdown  that  occurs  during  the 
leading  edge  of  a  switching  impulse  may  be  a  suitable  substitution  for  this 
type  of  measurement,  however. 

6.2  RATE  OF  VOLTAGE  CHANGE  AT  THE  END  OF  A  VERTICAL  HIGH-RESISTANCE 
CONDUCTING  CHANNEL  CREATED  BY  THE  IGNITION  OF  A  SUSPENDED 
DETONATING  FUSE 


It  has  been  suggested  that  lightning  could  be  triggered  by  suspending 
a  long  section  of. grounded  detonating  fuse  from  a  balloon  and  igniting  it 
(Ref.  16).  The  channel  created  by  the  ignition  has  a  resistance  of  about 
2  kft/m  for  a  duration  of  7  ms.  The  question  addressed  here  is  whether  or 
not  this  is  sufficient  time  for  the  upper  end  of  the  fuse  to  reach  a  poten¬ 
tial  that  is  nearly  the  same  as  the  ground  potential. 


The  equation  for  the  potential  of  a  transmission  line  is 

ii-  Lci^  +RC^  (37) 

3Z2  3t2 

where  L,  R,  and  C  are  assumed  to  be  slowly  varying  in  time  and  space.  If 
R  >>  ojL  for  frequencies  uj  of  interest,  the  first  term  on  the  right-hand 
side  of  Equation  37  may  be  neglected. 

=  RC  It  (38) 

az2  at 

Equation  38,  the  one-dimensional  diffusion  equation,  may  be  solved  using 
the  technique  of  separation  of  variables.  The  length  of  the  conducting 
cnannel  i  is  small  compared  to  the  distance  separating  the  earth  and  the 
charge  pocket  in  the  cloud;  therefore,  C  is  approximately  constant  over  the 
length  of  the  channel  (Ref.  17). 


The  initial  condition  is  V0(z)  3  Uz/Z  ana  the  boundary  conditions 
are  V(U,  t)  =  0  and  V(Z,  t)  =  0.  The  steady-state  solution  is 
V(z,  t  =  »)  =  0.  Separation  of  variables  takes  the  form  V  =  p(t)  V0 ( z ) , 


(39) 


Equations  for  V0  and  p  are 


dp  n2p 
dt  "  K 


(40) 


where  n  is  the  separation  constant.  By  expanding  V0  as  a  Fourier  series 
and  solving  Equation  40,  a  series  solution  is  obtained. 


*  V  *  2a  l  — jp-  sin  (rnrz/t)  exp  -(n2!^/^) 
n=l 


(41) 


where  a  =  u/ir.  The  time  constant  for  the  exponential  decay  of  the  end  of 
the  channel  to  zero  potential  is  z2RC/k2  *  (500  m) 2  (2xl03n/m)x 
(5xl0*12f/m)/ir2  =  0.25  ms.  Therefore,  the  lowest  Fourier  component  will 
decay  to  e"l  of  its  initial  amplitude  in  0.25  ms.  Higher  Fourier  com¬ 
ponents  will  decay  much  faster.  Therefore,  there  was  plenty  of  time  for 
the  upper  end  of  the  channel  to  reach  essentially  zero  potential,  and  the 
charging  time  required  does  not  provide  an  explanation  of  the  inability  to 
trigger  lightning  in  this  manner.  The  current  in  the  channel  is 


I  s  -  i  4y  =  2D  ^  (-I)n  COS(mr  z/  Z)  exp-(  n  2nZt/  Z2RC ) 

R  n=l 


142) 


where  0  =  itt/zk.  Equation  42  can  be  summed  analytically, 

I  =  O[04(u,q)-1]  (43) 

where  94  is  the  theta  function,  u  *  nz/2i,  and  q  =  exp-(  n^/ i2RC). 


The  time  interval  for  which  the  detenoting  fuse  retained  its  conduc¬ 
tivity  was  adequate  for  the  potential  of  its  upper  end  to  be  reduced  nearly 
to  ground  potential;  therefore,  the  failure  of  the  detonating  fuse  to  trig¬ 
ger  ligntning  can  not  oe  explained  by  an  inadequate  rate  of  charge  trans¬ 
fer  . 


Although  this  analysis  shows  that  the  resistance  of  the  channel 
formed  by  the  detonating  cord  was  low  enough  to  produce  electric  fields 
adequate  for  the  formation  of  an  upward  leader,  additional  analysis  is 
necessary  to  determine  if  the  resistance  was  too  high  to  sustain  the  leader 
growth.  This  is  a  possible  explanation  of  the  failure  to  trigger  light¬ 
ning.  It  is  also  possible  that  the  conducting  channel  was  sufficient  but 
the  atmospheric  conditions  were  not  suitable  for  lightning  formation. 
Lightning  is  not  triggered  by  every  wire-trailing  rocket  that  is  fired  into 
a  thunderstorm  at  the  North  Baldy  site. 
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